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Responsive  space-lift  is  a  new  trend  in  space  launches.  It  requires  the  guidance  system  of  the  launch 
vehicle  to  be  easily  prepared,  precise,  robust,  and  adaptive.  The  stellar/inertial  integrated  guidance 
system,  which  combines  two  distinct  types  of  navigation  data,  has  been  successfully  used  for  submarine- 
launched  ballistic  missiles.  This  paper  focuses  on  its  application  for  responsive  launch  vehicles.  Only  the 
prelaunch  positioning  and  orientation  errors  are  considered  in  order  to  provide  a  simple  illustration  of 
the  basic  principles,  which  can  be  conveniently  generalized  to  practical  launch  activities  with  a  well- 
developed  optimum  correction  method.  Error  transition  and  measurement  equations  are  constructed 
and  two  injection  error  correction  schemes  are  proposed.  Because  the  integrated  guidance  accuracy 
varies  with  regard  to  the  stellar  direction,  an  analytical  method  is  addressed  to  determine  the  optimal 
navigation  star.  Finally,  numerical  studies  are  conducted  and  the  results  verify  the  integrated  guidance 
policy.  Stellar/inertial  integrated  guidance  can  simplify  prelaunch  preparation  and  improve  the  orbital 
injection  accuracy.  Thus,  it  provides  an  alternative  that  does  not  rely  on  man-made  signals  for  responsive 
launch  vehicles. 
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1.  Introduction 

Over  half  a  century  of  spaceflight,  the  activity  of  entering  space 
is  no  longer  astonishing.  However,  a  reliable  space  launch  requires 
exhaustive  prelaunch  preparation,  which  is  time-consuming  and 
costly  work.  Consequently,  responsive  space-lift  has  become  a  new 
trend  in  spaceflight  [23,16].  The  objectives  of  responsive  space-lift 
include  low  launch  cost,  easy  preparation,  and  convenient  oper¬ 
ation,  which  could  make  space  launches  a  daily  activity.  Thus, 
guidance  systems  of  launch  vehicles  need  to  be  easily  prepared, 
precise,  robust,  and  adaptive. 

Many  guidance  approaches  have  been  investigated  since  the 
1960s,  all  of  which  need  to  measure  the  real-time  motion  states 
of  the  vehicle  to  validate  convergence  and  accuracy  [10,8,25,12], 
The  most  common  implementation  is  the  inertial  navigation  sys¬ 
tem  (INS).  However,  inertial  components  (the  gyroscopes  and  ac¬ 
celerometers)  used  in  an  INS  have  some  unavoidable  errors  with 
stochastic  properties.  For  vehicles  with  short  flight  times,  such  er¬ 
rors  might  be  acceptable.  For  long-time  duration  missions,  these 
errors  tend  to  increase  over  time  with  unbounded  position  error 
growth.  Thus,  for  precision  guidance,  the  INS  has  to  be  augmented 
with  external  aids,  among  which  the  most  extensively  and  suc¬ 
cessfully  used  is  the  Global  Positioning  System  (GPS)  [17,4,13],  Al¬ 
though  GPS  can  offer  highly  accurate  position  information  at  very 
low  cost,  the  likelihood  that  GPS  may  be  denied  to  the  user  ow¬ 
ing  to  interference  or  intentional  jamming  is  of  great  concern  and 
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thus  motivates  the  investigation  of  alternatives.  Typical  aid  devices 
include  a  vision  sensor  [24],  magnetometer  [11],  gravimeter  [19], 
and  star  sensor  [20,1,18,7],  In  the  case  of  employing  a  star  sen¬ 
sor,  a  Kalman  filter  has  been  used  to  combine  stellar  observations 
with  INS  outputs  to  enhance  the  navigation  performance  [20,1,18, 
7],  This  paper  proposes  another  data  fusion  method  for  the  ascend¬ 
ing  guidance  in  which  the  star-observing  window  is  so  short  that 
no  filter  can  be  used. 

The  stellar/inertial  integrated  system  is  a  simple  and  effec¬ 
tive  guidance  implementation  which  has  been  successfully  used 
for  submarine-launched  ballistic  missiles.  The  stellar-updated  iner¬ 
tial  guidance  system  incorporates  gyroscopes,  accelerometers,  and 
a  star  sensor  in  an  inertially  fixed  stable  element  that  mounted 
in  a  four-axis  or  five-axis  servo-stabilized  gimbal  system  (called 
inertial  platform  herein)  [14],  Prior  to  launch,  the  star  sensor  is 
pointed  in  a  computed  direction  through  a  leveling  and  alignment 
procedure,  and  the  attitude  is  maintained  during  launch.  At  the 
end  of  the  boost  stage,  prior  to  final  cutoff,  the  star  sensor  carries 
out  star  observation.  The  output  is  used  to  compute  the  target  im¬ 
pact  error  of  the  trajectory,  and  an  appropriate  velocity  increment 
is  made  to  reduce  the  predicted  error  to  zero  [15,6].  Thus,  the 
stellar/inertial  guidance  system  can  reduce  dependence  on  pre¬ 
launch  preparation,  enhance  performance  of  inertial  navigation, 
and  improve  guidance  accuracy,  all  of  which  are  in  accordance 
with  the  requirements  of  responsive  space-lift  [23],  Although  it 
has  not  been  adopted  by  traditional  launch  vehicles  mainly  owing 
to  a  lower  guidance  accuracy  requirement  compared  with  that  of 
ballistic  missiles,  it  is  a  suitable  guidance  solution  for  responsive 
launch  vehicles. 
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The  fundamental  of  the  guidance  system  is  that  observed  stel¬ 
lar  vector  is  used  to  determine  misalignment  angles  of  the  inertial 
platform,  from  which  guidance  errors  are  estimated  and  corrected. 
Therefore,  the  procedure  is  basically  a  problem  of  attitude  determi¬ 
nation  through  vector  observation.  This  kind  of  problem  was  first 
addressed  and  investigated  by  Wahba  in  the  1960s  [21],  and  many 
algorithms  have  been  developed  since  then  [2,9,5],  Most  require 
two  observed  vectors  to  achieve  convergence,  but  some  function 
well  with  a  long-term  observation  of  one  vector  [22],  Similarly,  ac¬ 
cording  to  the  number  of  vectors  being  measured,  stellar/inertial 
integrated  guidance  systems  can  be  categorized  as  single-star  and 
double-star  schemes.  The  difference  is  that  impact  point  down- 
range  and  cross-range  deviations  can  directly  relate  to  the  outputs 
of  star  sensor  and  there  is  no  need  to  determine  the  attitude  of 
the  inertial  platform.  For  there  are  only  two  accuracy  indices,  the 
single-star  scheme  can  theoretically  achieve  the  same  accuracy  as 
the  double-star  scheme  by  observing  a  specific  vector,  which  is 
called  the  optimal  stellar  direction.  For  the  launch  vehicle  guid¬ 
ance  system,  there  are  more  than  two  accuracy  indices,  and  that 
results  in  differences  between  the  two  schemes.  In  spite  of  that, 
the  structural  configuration  of  gimbals  system  and  the  correction 
procedure  of  orbital  injection  errors  in  the  single-star  scheme  are 
simpler  and  more  feasible,  and  thus,  this  paper  explores  the  fun¬ 
damentals  and  applications  of  the  single-star  scheme  rather  than 
those  of  the  double-star  scheme. 

During  a  practical  space  launch,  misalignment  angles  are  in¬ 
troduced  by  many  uncertainties,  including  prelaunch  positioning 
and  orientation  errors,  initial  alignment  and  leveling  errors,  and 
inertia  component  errors.  The  stellar  vector  measurement  cannot 
distinguish  these  errors  into  classes,  and  thus,  a  synthetic  correc¬ 
tion  based  on  optimum  estimation  theory  is  required,  supposing 
that  the  error  distributions  are  known  [15,6],  Because  the  opti¬ 
mum  correction  method  has  been  well  developed,  for  the  sake  of 
simple  illustration,  only  prelaunch  positioning  and  orientation  er¬ 
rors  are  taken  into  account  in  this  paper.  It  is  also  assumed  that 
a  five-axis  servo-stabilized  gimbal  system  is  adopted  so  that  the 
optical  axis  of  the  star  sensor  can  aim  at  an  arbitrary  navigation 
star  [6],  Undoubtedly,  a  strap-down  star  sensor  has  some  virtues 
over  the  gimbal  system.  We  adopt  the  servo-stabilized  star  tracker 
mainly  for  its  maturity  in  application,  and  the  principle  of  the  stel¬ 
lar/inertial  integrated  guidance  using  a  strap-down  star  sensor  is 
almost  the  same  with  that  studied  in  our  paper. 

The  paper  is  organized  as  follows.  In  Section  2,  error  transition 
and  measurement  equations  are  derived  and  two  orbital  injection 
error  correction  schemes  are  proposed.  An  analytical  method  to 
determine  the  optimal  stellar  direction  is  presented  in  Section  3, 
and  in  Section  4,  numerical  studies  are  conducted  to  validate  the 
foregoing  conclusion. 

2.  Integrated  guidance  modeling 

2.1.  Error  transition  equation 


Prelaunch  positioning  errors  can  be  represented  by  orthogonal 
coordinate  errors  Axo,  Ayo.  and  Azo  in  the  geocentric  fixed  frame, 
or  the  longitude  error  Aao,  geographic  latitude  error  ABo,  and  al¬ 
titude  error  AHo  [3],  Using  spherical  triangle  equations,  we  easily 
obtain  the  relationship  between  the  two  groups  of  errors 
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where  A0,  Bo,  N,  and  H0  are  the  launch  azimuth,  geographic  lat¬ 
itude,  radius  of  the  prime  vertical,  and  altitude,  respectively.  Be¬ 
cause  the  altitude  error  does  not  affect  the  inertial  platform  align¬ 
ment  and  leveling,  it  cannot  be  modified  through  stellar  measure¬ 
ment.  Consequently,  only  AAo,  ABo,  and  AAo  (the  launch  azimuth 
error,  representing  the  prelaunch  orientation  error)  are  considered 
in  the  following  discussion. 

Deviations  of  the  motion  states  induced  by  the  positioning  and 
orientation  errors  can  be  divided  into  two  parts:  the  initial  er¬ 
rors  of  navigation  computation  and  the  integral  errors  brought  by 
the  misalignment  between  the  inertial  platform  frame  P  and  the 
launch  point  inertial  frame  I.  The  axes  of  the  coordinate  system 
P  coincide  with  the  gimbal  axes  [6],  The  x  axis  of  frame  I  is  ori¬ 
ented  along  the  range,  and  the  z  axis  is  perpendicular  to  the  local 
horizontal  plane.  For  the  former  part  of  the  deviations,  the  initial 
position  error  is 
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and  the  initial  velocity  error  is 


Av(t0)  =  ft>ex  •  Ar(t0),  (3) 

where  (oex  is  the  anti-symmetric  matrix  of  the  Earth  spin  angular 
velocity  toe.  The  definition  of  toex  is 
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The  integral  errors  can  be  described  in  terms  of  the  three  mis¬ 
alignment  angles  ax,  ay,  and  otz.  Suppose  that  frame  I  successively 
rotates  about  the  three  axes  with  angles  ax,  ay,  and  olz  to  coincide 
with  frame  P.  With  first-order  approximation,  the  direction  cosine 
matrix  for  transformation  between  the  two  frames  is  expressed  as 

I-  1  -«z  *y  1 
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The  direction  cosine  matrix  for  transformation  between  the  geo¬ 
centric  inertial  frame  E  and  frame  I  is 


C'E  =  M2[-(90°  +  A0)]  •  Mi [B0]  •  M3[-(90°  -  A0)],  (6) 

where  JW,  (i  =  1,2,3)  are  elementary  rotation  matrices  [3],  Fur¬ 
thermore,  because  the  inertial  platform  is  aligned  and  leveled  at 
the  actual  launch  point,  taking  no  account  of  the  alignment  and 
leveling  errors,  we  obtain  the  direction  cosine  matrix  for  transfor¬ 
mation  between  frame  P  and  frame  E 


CPE  =  M2[-(90°  +  Ao  +  AA0)]  •  Mi [B0  +  AB0] 

x  M3[— (90°  —  Ao  —  AAo)].  (7) 

Similarly,  with  first-order  approximation,  we  can  obtain  another 
expression  for  the  matrix  C{,  by  multiplying  Eqs.  (6)  and  (7) 
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Comparing  Eq.  (5)  with  Eq.  (8),  we  find  the  relations  between  the 
misalignment  angles  and  the  positioning  and  orientation  errors  are 

[ax  "1  To  —  sin  Ao  cos  Ao  cos  Bo  1  ["  AAo"l 

ay  =  -1  0  sinBo  AB0 

OLz  J  |_  0  -cos Ac  -cosBosinAo  J  |_  AAo  J 

=  P  A. 


Since  the  optical  axis  OXB  is  nearly  collinear  with  the  stellar  direc¬ 
tion  vector,  it  is  reasonable  to  assume  that  the  cosine  value  of  the 
included  angle  is  1.  If  the  outputs  of  the  star  sensor  are  $  and  rj, 
then  vector  S  in  frame  B  can  be  expressed  as 

SB=[1  -$  -r/]T.  (13) 

(9)  Therefore,  the  expression  of  vector  S  in  frame  P  is 


The  matrix  P  is  invertible  when  B  /  90°;  hence,  generally  a 
uniquely  corresponds  to  A.  In  other  words,  A  can  be  uniquely  de¬ 
termined  by  measuring  a. 

According  to  the  derivation  above,  the  relationship  between  the 
integral  errors  and  A  can  be  derived,  although  the  expression  is 
rather  complicated.  If  we  define  Rk  and  Vk  to  be  the  vehicle  po¬ 
sition  and  velocity  vectors  when  the  star  alignment  is  carried  out, 
then  expressions  of  the  motion  state  deviations  in  terms  of  the  po¬ 
sitioning  and  orientation  errors  can  be  derived.  They  are  expressed 
in  a  compact  form  as 

r  ARk  1  _  r  dRk/dAo  8Rk/8B0  8Rk/8X0  1  [ 
[AVk\-[dVk/8A0  8Vk/8B0  8Vk/8X0\ 

=  C  •  A.  (10) 

The  matrix  C  can  be  computed  in  advance  for  that  particular  po¬ 
sition. 


Sp  =  CpbSb.  (14) 

According  to  the  transfer  relationship  between  coordinate  sys¬ 
tems,  we  obtain  another  expression  of  vector  S  in  frame  B 

Sb  =  Cb-Cp -S,.  (15) 


Inserting  the  expressions  of  the  related  matrices  and  vectors  into 
the  preceding  equation,  the  measurement  equation  is  obtained  as 


Vector  w  =  [f ,  rj]T  is  adopted  to  denote  the  outputs  of  the  star 
sensor  in  the  following. 


2.2.  Measurement  equation 

We  denote  the  elevation  angle  and  azimuth  angle  of  the  stel¬ 
lar  direction  unit  vector  S  in  the  frame  I  by  es  and  crs,  as  shown 
in  Fig.  1.  The  elevation  angle  is  measured  upward,  while  the  az¬ 
imuth  angle  is  measured  clockwise.  Thus,  vector  S  in  frame  I  can 
be  described  as 

Si  =  [  cosescosors  sines  cosessinors  ]r.  (11) 

To  point  to  the  navigation  star,  the  star  sensor  body  frame  B 
rotates  with  regard  to  frame  P,  under  the  assumption  that  frame 
P  coincides  with  frame  I  after  alignment  and  leveling.  Thus,  frame 
P  firstly  rotates  about  the  y  axis  by  —  as  and  then  about  the  z  axis 
by  es  to  coincide  with  frame  B,  as  shown  in  Fig.  2.  The  direction 
cosine  matrix  for  the  two  frames  is  then  derived  as 

CBP  =  M3(es)M2(-as) 

cos  es  cos  <rs  sines  cosessincrs  "1 
—  sin  es  cos  crs  coses  —  sinessin<Ts  .  (12) 
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2.3.  Correction  of  orbit-injection  error 

When  launching  a  satellite,  six  orbital  element  constraints  are 
usually  specified  in  order  to  accomplish  scheduled  tasks.  The  ex¬ 
pected  right  ascension  of  ascending  node  is  generally  obtained  by 
choosing  launch  window,  while  the  true  anomaly  is  obtained  by 
in-orbit  phasing.  Constraints  on  the  other  four  orbital  elements, 
namely  semi-major  axis  a,  eccentricity  e,  orbital  inclination  i,  and 
argument  of  perigee  to,  usually  hold  by  iterating  launch  program. 
Sometimes  another  four  parameters,  usually  the  geocentric  dis¬ 
tance  rk,  magnitude  of  velocity  vk,  flight  path  angle  &k,  and  orbital 
inclination  i,  are  used  as  the  iterating  variables  instead.  The  single¬ 
star  scheme  has  only  two  measurements.  Although  some  newly 
developed  filtering  algorithms  can  get  the  three  attitude  angles, 
the  procedure  of  star  alignment  takes  only  a  few  minutes  within 
which  no  filter  could  converge.  Thus,  to  achieve  optimal  accu¬ 
racy,  two  of  the  four  iterating  variables  should  be  selected.  The 
double-star  scheme  is  capable  of  determining  the  three  misalign¬ 
ment  angles  of  the  inertial  platform,  and  hence  all  four  parameters 
can  be  corrected,  resulting  in  more  promising  accuracy.  However, 
noticing  that  only  three  control  variables  (direction  of  thrust  vec¬ 
tor  and  the  time  of  cutoff)  are  available  during  the  correction,  if 
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four  injection  parameters  are  corrected  simultaneously,  an  opti¬ 
mal  control  algorithm  should  be  employed  to  reach  a  satisfying 
result.  For  the  correction  duration  is  not  long,  simulations  show 
that  the  convergence  of  the  guidance  algorithm  is  unstable.  If  only 
two  parameters  are  corrected,  the  known  and  simple  velocity-to- 
gain  guidance  algorithm  can  be  used  with  reliable  convergence. 
Furthermore,  the  structural  configuration  of  a  gimbal  system  con¬ 
taining  two  star  sensors  is  much  more  complicated.  Therefore,  the 
single-star  scheme  is  recommended  and  investigated  in  this  pa¬ 
per. 

A  number  of  combinations  of  two  members  of  the  iterating 
variables  can  be  chosen  as  correction  parameters,  among  which 
two  schemes  are  suggested  in  consideration  of  the  importance  of 
injection  requirements.  One  is  the  semi-major  axis  a  and  orbital 
inclination  i,  and  the  other  is  the  flight  path  angle  &k  and  or¬ 
bital  inclination  i.  It  is  known  that  a  is  a  scalar  that  represents  the 
orbital  energy,  while  the  measurements  of  the  star  sensor  indi¬ 
cate  the  inertial  platform  drift,  which  somewhat  is  a  vector.  Thus, 
we  deduce  that  the  integrated  system  may  not  obviously  improve 
the  accuracy  of  a.  Contrarily,  &k  and  i  both  imply  information  of 
the  orientation,  and  in  a  manner  can  be  respectively  analogized  to 
the  longitudinal  range  and  lateral  range  of  ballistic  missiles.  Thus, 
a  distinct  correction  effect  can  be  expected.  Although  the  correc¬ 
tion  difference  between  the  parameters  vanishes  when  only  the 
positioning  and  orientation  errors  are  taken  into  account,  the  de¬ 
duction  is  verified  by  numerical  simulations  (results  not  presented 
here)  when  all  the  related  errors  that  affect  the  misalignment  an¬ 
gles  are  involved. 

Partial  derivatives  of  the  correction  parameters  with  respect  to 
the  motion  states  are  derived  as  follows.  We  first  define  the  direc¬ 
tion  cosine  matrix  for  transformation  between  the  frame  /  and  the 
launch  point  equatorial  inertial  frame  El  (with  the  x  axis  point¬ 
ing  toward  the  intersection  of  the  equator  and  the  meridian  of  the 
launch  point)  as 


Similarly,  by  taking  derivatives  of  Eqs.  (21)  and  (22),  the  partial 
derivatives  of  i  and  0k  with  respect  to  the  cutoff  states  are  derived 
as 
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and 
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We  denote  the  correction  parameters  of  the  two  schemes  pro¬ 
posed  above  as  m i  and  m2.  Combining  Eq.  (18)  with  Eqs.  (23)-(25), 
we  obtain  the  partial  derivatives  expressed  in  frame  I 
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cf 1  [-B0]  •  M2  +  A0  j .  (17) 

If  we  denote  the  geocentric  vector  of  the  launch  point  by  rLl  it  is 
clear  that  the  position  and  velocity  vectors  before  star  alignment 
in  frame  El  are 

J*a  =  cf‘(*  +  rO.  (18) 

l  1/Efc  =  cf^  vfc. 

The  angular  momentum  is 


hk  =  REk  X  V  Ek. 


(19) 


According  to  the  relations  between  the  orbital  elements  and  mo¬ 
tion  states  [3],  one  has 


a  _  REk 


2-vk 

_1  hkz  , 
t  =  cos  — ,  and 
hk 


©k  =  sin' 


REk  '  V  Ek 
REkVEk 


(20) 

(21) 

(22) 


where  vk  —  and  /x  is  the  gravitational  constant  of  the 

Earth. 

Taking  the  derivative  of  Eq.  (20),  the  partial  derivative  of  a  with 
respect  to  the  cutoff  states  is  obtained  as 


da 

2/x2 

REk 

dREk 

(REkV2Ek-2^2 

REk  ’ 

da 

2  nR2EkvEk 

vEk 

3 V7k 

(REkV2k  -  2/x)2 

VEk 

From  Eqs.  (10)  and  (26),  the  relationship  between  the  motion  state 
deviations  and  the  initial  errors  is  obtained  as 


From  Eq.  (9),  we  obtain  A  =  P_1  ■  a  =  D  ■  a.  Inserting  the  expres¬ 
sion  into  the  foregoing  equation,  we  express  the  injection  errors  in 
terms  of  the  misalignment  angles  in  a  compact  form  as 


I"  Ami  1  _  I"  dm\/dux  dm^/doiy  3mi/3o!zl 
LAm2J-L3m2/3«x  9m2/ day  3m2/3a!zJ 

a  =  Na.  (28) 


3.  Optimal  navigation  star 

Eq.  (28)  indicates  that,  when  only  the  positioning  and  orienta¬ 
tion  errors  are  taken  into  account,  it  is  the  projection  as  of  the 
vector  formed  by  the  misalignment  angles  in  the  subspace  Ns 
spanned  by  ttj  and  n2  that  affects  the  injection  accuracy.  Hence, 
according  to  Eq.  (16),  if  the  subspace  Hs  —  span{h^h2\  is  designed 
to  be  equivalent  to  Ns  by  choosing  the  vectors  hi  and  h2  (basically 
the  direction  angles  es  and  as  of  the  navigation  star),  then  the  out¬ 
put  of  the  star  sensor  w  contains  all  useful  information.  The  star 
in  the  particular  direction  is  the  optimal  navigation  star. 
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Table  1 

Guidance  accuracy  comparison. 


Aa  /m  Mean 

MSE 

Ai  /deg  Mean 

MSE 

A©  /deg  Mean 

MSE 


Analysis  covariance 

Inertial  Integrated 

0.0  0.0 

1981.861  74.083 

0.0  0.0 

1.067  xlO-1  3.847  x  10-2 

0.0  0.0 

1.272  xlO-2  5.476  x  IQ-3 


Monte  Carlo  sampling 


-3.857  x  10“3 
1.092  x  liH 
-3.607  x  10  4 
1.296  x  10-2 


Integrated 

-2.684 

73.558 

-1.394  xl0“3 
3.819  x  10-2 
—1.984  x  10~4 
5.438  x  10“3 


We  can  derive  as  from  Eq.  (16) 


a s  =  Ht(HHt)  ’w.  (29) 

Consequently,  the  correction  equation  is 


The  corrected  injection  parameters  are  Ami  -  Amis  and  Am2  - 
Arri2s.  Because  the  orthogonal  complement  spaces  of  Ns  and  Hs 
should  be  equivalent,  the  optimal  stellar  direction  can  be  obtained 
by  solving 


h  i 


Ml  X 

|ni  xn2|' 


(31) 


4.  Numerical  simulation  analysis 


In  this  section,  numerical  simulations  conducted  to  evaluate 
the  effect  of  the  integrated  guidance  technique  are  presented.  It 
is  supposed  that  a  Sun-synchronous  satellite  with  an  orbital  al¬ 
titude  of  485  km  is  launched.  The  geographic  coordinates  of  the 
launch  point  were  set  to  111.61°E,  38.85°N,  and  the  launch  az¬ 
imuth  A0  was  192.94°.  The  positioning  and  orientation  errors  were 
assumed  to  have  normal  distributions  with  zero  means,  and  the 
mean  square  errors  (MSEs)  were  set  as  cta  =  [40"40"60"]r. 

After  some  skipped  calculations,  the  partial  derivative  matrix 
of  the  correction  parameters  with  regard  to  the  positioning  and 
orientation  errors  was  obtained 
'  3a/3 Ar 
3i/3AT 
_  3@fc/3Ar  _ 


'  5.195  x  106 
4.083  x  10-1 
1.587 


8.720  xlO6  7.933  xlO5' 
3.646  xlO"1  6.700  x  lfT1 
1.621  6.220 


A  star  in  the  direction  of  es  =  45°  and  crs  =  45°  was  adopted 
as  the  navigation  star.  The  analysis  covariance  method  and  Monte 
Carlo  sampling  were  employed  to  obtain  the  guidance  accuracy 
for  both  the  integrated  and  purely  inertial  guidance  schemes.  The 
number  of  samples  was  400.  Table  1  lists  the  results. 

The  results  in  Table  1  show  the  accordance  of  the  results  ob¬ 
tained  using  the  two  methods.  By  comparing  the  accuracy  data 
of  the  two  schemes,  it  is  clear  that  both  the  correction  schemes 
proposed  above  effectively  improve  the  injection  accuracy.  When 
compared  the  results  with  the  INS/GPS  system,  the  stellar/inertial 
guidance  system  is  less  accurate.  Navigation  data  compensation  is 
only  provided  after  cutoff,  while  the  INS/GPS  is  available  for  all 
the  boost  phase.  However,  our  intention  is  to  propose  an  alter¬ 
native  that  does  not  rely  on  any  man-made  signals  and  that  can 
significantly  improve  the  injection  accuracy. 

Figs.  3,  4,  and  5  show  the  integrated  guidance  accuracies  of  the 
semi-major  axis,  orbital  inclination,  and  flight  path  angle  with  re¬ 
gard  to  the  direction  of  the  navigation  star,  respectively.  The  results 


reveal  that  the  direction  of  the  navigation  star  strongly  affects  the 
guidance  accuracy. 

When  comparing  the  results  presented  in  the  figures  and  Ta¬ 
ble  1,  it  is  found  that  the  accuracy  of  integrated  guidance  is  some¬ 
times  even  poorer  than  that  of  the  purely  inertial  guidance.  The 
reason  is  that  the  correction  coefficients  used  in  the  simulations 
are  defined  in  Eq.  (30)  with  regard  to  the  optimal  navigation  star. 
When  the  adopted  star  considerably  departs  from  the  optimal  one, 
the  original  inertial  deviation  may  be  overly  corrected.  To  avoid 
this  problem,  different  correction  coefficients  corresponding  to  the 
adopted  stars  should  be  used  during  simulations. 

Fig.  6  consists  of  sections  of  Fig.  3.  The  plots  show  that  two 
stars  will  have  the  same  correction  effects  when  they  are  collinear 
in  terms  of  the  origin  of  frame  /;  i.e.,  esi  =  — eS2  and  |crsi  —  crs2|  = 

Employing  the  method  presented  in  Section  3,  the  optimal 
stellar  direction  was  obtained.  For  the  scheme  using  the  semi- 
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major  axis  and  the  orbital  inclination,  the  results  are  eSOpti  = 
60.920°  and  crsopti  =  84.795°.  For  the  second  scheme,  the  results 
are  eSOpt2  =  37.393°  and  orSopt2  =  74.544°.  Computations  show 
that  the  original  inertial  injection  errors  can  be  totally  corrected 
through  measurements  of  the  optimal  navigation  star.  Moreover, 
as  mentioned  above,  the  use  of  another  star  that  is  collinear  with 
the  optimal  star  in  terms  of  the  origin  of  frame  /  can  achieve  the 
same  accuracy. 

Employing  the  optimum  correction  method,  more  simulations 
were  conducted  accounting  for  all  related  errors.  The  results  show 
that,  when  the  optimal  navigation  star  is  adopted,  the  integrated 
guidance  can  remarkably  improve  the  injection  accuracy.  If  the 
semi-major  axis  and  the  inclination  are  corrected,  the  accuracies 
are  approximately  enhanced  by  1.5  times  and  5  times.  In  the 
second  correction  scheme,  the  accuracies  of  the  inclination  and 
flight  path  angle  are  approximately  enhanced  4  times  and  4.5 
times. 

5.  Conclusions 

Referring  to  the  principles  of  stellar/inertial  integrated  guidance 
used  for  ballistic  missiles,  this  paper  studied  the  application  of 
such  guidance  to  responsive  launch  vehicles.  Error  transition  and 
measurement  equations  were  derived  when  only  the  positioning 
and  orientation  errors  were  taken  into  account.  Two  injection  error 
correction  schemes  were  proposed,  and  numerical  analyses  were 
conducted  to  evaluate  the  performance  of  the  integrated  guidance 
technique.  Results  show  that  the  stellar/inertial  guidance  system  is 
a  suitable  alternative  for  responsive  space-lift. 

Nevertheless,  only  the  fundamentals  of  the  approach  are  dis¬ 
cussed  here.  The  elements  that  affect  the  realization  of  an  inte¬ 
grated  guidance  system  are  simplified,  and  the  correcting  proce¬ 
dure  has  not  yet  been  implemented.  Thus,  many  issues  still  require 
investigation. 
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